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THREE ANALOGUES OF STERN’S DIATOMIC SEQUENCE 


SAM NORTHSHIELD 


Abstract. We present three analogues of Stern’s diatomic sequence. 


1. Introduction 

Stern’s diatomic sequence a\ = 1, a^n = « 2 n+i = a n + a n+ \ is a particularly well studied 

sequence (see, e.g., in. 0, m and references therein, as well as m )■ The first section 
is devoted to showing that this sequence is interesting. In particular, we shall look at the 
following properties. 


• n i-A a n+ i/a n is a bijection between the positive natural numbers and the positive 
rational numbers, 

• n/2 fc i—>■ a n /a n+ 2 k extends to a continuous strictly increasing function on [0,1] known 
as “Conway’s box function” (it’s inverse is ?(x), Minkowski’s question-mark function), 

• It shares a number of similarities to the Fibonacci sequence; in particular, it has a 
Binet type formula. 


The remaining three sections are devoted to three analogues of Stern’s sequence: 

• We replace addition by another binary operation; in particular, we define b\ = 0, b 2 n = 
b n , & 2 n+i = b n © b n .|_i where x © y = x + y + \/^xy + 1. This sequence is related 
to Stern’s sequence and arises from certain sphere packings. It has apparently not 
appeared before in the literature. 

• We replace the complementary indexing sequences {2n} and {2n + 1} by another pair 
of complementary sequences; in particular, let R\ = 1, R a ( n ) = Rn-,Rp(n) = Rn + R n +1 
where a(n) = [n4> — l/</> 2 _|, /3(ra) = [ncf> 2 + J form a specific pair of complementary 
Beatty sequences. This sequence has been extensively studied as R. n is the number of 
ways n can be represented as a sum of distinct Fibonacci numbers. 

• The known Binet type formula for Stern’s sequence [9] is written in terms of the 
sequence S 2 (n) (:= the number of terms in the binary expansion of n). We replace 
S 2 (n)by Si7’(n)(:= the number of terms in the Zeckendorf representation of n). This new 
sequence, apparently not studied before, is an integer sequence with several interesting 
(and several conjectural) properties. 
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2. Stern’s Diatomic Sequence 


Consider the following “diatomic array” [1] formed as a variant of Pascal’s triangle; each 
entry is either the value directly above or else the sum of the two above it. 

1 1 
1 2 1 
1 3 2 3 1 

1 4 3 5 2 5 3 4 1 

15473857275837451 


The word “diatomic” is used here since every entry of the diatomic array gets its value from 
either one or two entries above and gives that value to three entries below, hence has “valence” 
4 or 5 (hence the diatomic array models a kind of crystalline alloy of two elements). 

Ignoring the right most column and reading the numbers as in a book, we get Stern’s 
diatomic sequence: 

1,1,2,1,3,2,3,1,4,3,5,2,5,3,4,1,5,... 

The sequence is thus defined by the recurrence 

®1 — lj 0'2n — On> ®2n+l — T Un+1- (1) 

We define ao to be 0 (the value consistent with gq.o+i = ao + «i)- 

Perhaps the most celebrated property of this sequence is that every positive rational number 
is represented exactly once as a n+ \/a n . See, for example, (3j or [9]. We rephrase this fact as 
a theorem. 

Theorem 2.1. Every ordered pair of relatively prime positive integers appears exactly once in 
the sequence {(a n ,a n +i)}. 

Proof. For an ordered pair, consider the process of subtracting the smallest from largest (stop 
if equal). For example, (4, 5) eA (4,1) and (7, 3) e-)- (4, 3). By the definition of Stern’s sequence, 

(fl2n j ®2n+l) j (®2n+l j ®2n+2) 1 t (fl n , Un+l). 

Every relatively prime pair appears (if not, then there is an ordered pair not on the list with 
lowest sum. Apply the process; the result has lower sum and so is (a n , a n+ 1 ) for some n and so 
the original pair is either ((Qn, « 2 n+i) or (®2n+i> R 2 n+ 2 ))- Every relatively prime pair appears 
exactly once since, if not, then there exist m < n with (a m ,a m +i) = (a n ,a n _|_i) and such that 
m is as small as possible. Applying the process to both implies \_m/2\ = \n/2\ and thus 
a m = a m+ 1 = o m+ 2 which is impossible. □ 

One can then rewrite any sum over relatively prime pairs in terms of Stern’s sequence. As 
an example, we rephrase the Riemann hypothesis. First note that n i— > GQn/cQn+i is an 
explicit bijection from Z + to Q 0 (0,1). Then the Riemann hypothesis is equivalent to 

2 iria 2 n/ci‘ 2 n +1 = 0(x 1/2+e ). 

^2n+l 

Briefly why this is so: the Mobius function can be written as p(ji) := Yli<k<n gcd(?i fc)=i e 2mk / n 
and so the left hand side is really just Merten’s function M(x) := Yh n <x / i ( n )- The connection 
between Merten’s function and the Riemann hypothesis is well-known; see for example [5]. 
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Minkowski’s question mark function was introduced in 1904 as an example of a “singular 
function” (it is strictly increasing yet its derivative exists and equals 0 almost everywhere). It 
is defined in terms of continued fractions: 


?(*) = 2]T 

n= 1 


(_l)n+i 

2&1 +a2 + ---+fln 


where x = l/(oi + l/(a 2 + 1/(03 + ...)))■ By Lagrange’s theorem that states that the continued 
fraction representation of a quadratic surd must eventually repeat, it is clear that ?(x) takes 
quadratic surds to rational numbers. 

The function 



&k 

a 2 n +k 


extends to a continuous strictly increasing function on [0,1]. This function is known as “Con¬ 
way’s box function” and its inverse is Minkowski’s question mark function ?(x). See [9] for a 



Figure 1. The graphs of y = f(x) and its inverse y =l{x). 


proof. 

The functions f(x) and ?(x) extends to homeomorphisms (or, equivalently, are restrictions 
of homeomorphisms) between two fractals. 




Figure 2. Sierpinski gasket and an Apollonian circle packing 
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Stern’s sequence is related to the Fibonacci sequence in a number of ways. For example, 


the Fibonacci sequence is 

a subsequence: 



m 



1 

i 


0 

1 

i 

0 

2 

3 1 

1 4 

3 GO 

2 5 

3 4 1 

1 5 4 

7 3 0 5 

7 2 7 5 

8 3 7 4 5 1 


It is easy to see that 

a J(n) = F n (2) 

where J(n) := ( 2 n — (—l) n )/3 is the Jacobsthal sequence [T2], A001405]. 

A new result by Coons and Tyler |3j identifies and proves the asymptotic upper bound: 

The constants involved in this formula are perhaps not so surprising since, by formula (2), it 
is clear that 

hm a J(n) /((3J(n)) 10 ^) = 1/^5. 

n—too v 7 

Stern’s sequence has a few remarkable similarities to the Fibonacci sequence (see [9] and 
[lOj). For example, Stern’s sequence satisfies a modified Fibonacci recurrence: 

0*71 +1 — T ®n—l 2 (Q, n —i mod CL n j. 

Next, certain diagonal sums across Pascal’s triangle yield the Fibonacci sequence while the 
corresponding sums across Pascal’s triangle modulo 2 yield Stern’s sequence: 


1 1 

1 2 0 

1 0 3 

[T] 4 6 4 


1 0 0 


1 

1 1 _ 

i o on 

i0ii 


0 0 


Fn+1 ~ JO 


2 i-\-j=n 


i+j 


Recall Binet’s formula 


a n +1 — 

2 i-\-j=n L 


i + j 


mod 2 


in+l 


—rn +1 n 


Fn +1 — 


4>~ 


k =0 


fc-rn—k 


Stern’s sequence satisfies a similar formula: 


a n+1 = JO (j s ^a s ^ n ~ k) 


k =0 


( 3 ) 


where a := (1 + — 3)/2 is a sixth root of unity and S 2 (n) is the number of ones in the binary 

expansion of n [12] A000120]: 


n 01234567 8... 
s 2 {n) 0 1 1 2 1 2 2 3 1... 
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Why this is true: If G(x) := 'f^a S2 ^ > x n and F(x) is the generating function for {a n+ 1 }, then 
G(x) = (1 + <jx)G(x 2 ) and F(x) = (1 + x + x 2 )F(x 2 ). For real x, since (1 + ax)(l + ax) = 
1 + x + x 2 , \G(x)\ 2 = F(x) and the result follows by equating coefficients. 


3. Replacing addition by another operation 
Definition 3.1. For non-negative real numbers a,b, let 

a ©6 = a + b + V 4a6 + 1 
a Qb = a + b — y/Aab + 1 

Proposition 3.2. If a, b, c, d > 0 and \ad — bc\ = 1 then (ac ) © ( bd ) = (a + b)(c + d). 

Proof. If (ad — be) 2 = 1, then 

(ad + be) 2 = 1 + 4 abed 

and thus 

(ac) © (bd) = ac + bd + V4 abed + 1 = ac + bd + ad + be. 

□ 

Remark 3.3. By the Fibonacci identity 

Fn-iF n+1 = F 2 + (-If, 

it follows that 

(F n —\F n ) © (F n F n+ 1 ) = (F n _i + F n )(F n + -Fn+i) = -F r) ,+il ? r) ,_|_ 2 . 

and so the sequence x n := F n F n+ i satisfies the modified Fibonacci recurrence 

%n +1 — © a: n _r. 


Here we define the first new sequence. 


Definition 3.4. Let b± 


0, and for n> 1, 

&2n = b n 

&2n+l = b n © &n+l- 


The sequence begins 


0,0,1, 0,2,1, 2, 0, 3, 2, 6 ,1, 6 , 2, 3, 0,4,3,10, 2,15,6,12,1,12, 6,15,.. 


It is not immediately clear that this sequence must always be integral. One way to show 
this is to express each bk as a product of elements of Stern’s sequence (Theorem 3.6, below). 
First we must prove a lemma. 


Lemma 3.5. For m,n>0, if m + n = 2 J — 1 then a m+ ia n+ i — a m a n = 1. 


Proof. We prove this by induction on j. If m+n = 1, then a m+ \a n+ \ — a m a n = 0102 — ao«i = 1 
and the result holds for j = 1. Suppose now that the result holds for a fixed j and that 
m + n = 2 J+1 — 1. Without loss of generality, m = 2k + 1 and n = 21 for some k, l > 0 (and 
so k + l = 2 J — 1). Then 

CLm+l a n+l — (Lma n = 02fc+2«2«+l — &2k+l a 2l 

= + O'2/u-i) — (a*, + afc + i)aj = afc + ia; + i — a^cq = 1 
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and the result follows. □ 

Theorem 3.6. If 2 J < k < 2- ?+1 , then 

bk = a 2 j+l_ k a k _ 2 j. 

Proof. If k = 2 J then, because ao = 0, bk = 0 = a 2 j+i_ k a k - 2 j = a 2 i-k a k- 2 i~ 1 - 

Let Xk := a 2 J+i-fc a fc_ 2 i where k E (2 J , 2 J+1 ). Then 2k, 2k + 1 E (2 J+1 ,2 J+2 ) and thus 


,T2fc — a 2 j+ 1 -2fe a 2fe-23 — a 2i-k a k-2i~ 1 — x k 

and, by lemma 3.5 and proposition 3.2, 

x 2fc+l = a 2J+ 1 -(2fe+l) a 2fe+l-2t 

= a2(23-fe-l)+l a 2(fe-2J-l)+l 
= («2J — fc—1 + a 2i-k) ' ( a fc-2J- 1 + a fe+l-2J- 1 ) 

= ( a 2J-fc-l a fc+l-2J- 1 ) © ( a 2J-fe a fc-2J- 1 ) = ^fc+l © 

Hence 6 ^ = for all fc, and the result follows. □ 

Corollary 3.7. b n € N. 

As seen in section 2 , Stern’s diatomic sequence leads to a construction of Conway’s box 
function f(x), the inverse of Minkowski’s question-mark function. The sequence {bk} gives 
rise to a similar function that turns out to be closely related to f(x). 

Definition 3.8. For k,n £ N, k < 2 n , let 

(k_\ , = _bk_ 

9 \2 n )' b 2 , +k ' 

Theorem 3.9. The function g[x) extends to a continuous function on [0,1] that satisfies, for 
x E (2 -J_1 , 2~ U, 

g(x) = /( 2 3+l x — 1)[1 — jf( 2x)\ 
where f(x) is Conway’s box function. 

Proof. Let x = k/2 n . Then 2 n ~i~ l < k < 2 n ~i for some j > 0. Since 2 n < 2 n + k < 2 n+1 , it 
follows from Theorem 3.6 that 


bk — a 2 n -i —k a k—2 n -i- 1 an d b 2 n +k — a 2 n_ k a k . 
By [9j formulas (2) and (3)], 

a 2 n -fc = JUk + «2 n -i-k 

and thus 

k \ _ b k _ a 2 n-j-k a k-2 n ~^- 1 

2" J _ ~ 


g{x) = g 


b 2 n +k a 2 n_ k a k 

( a 2 n -k — j a k) a k-2 n ~i~ 1 a k—2 n ~i~ 1 


a 2 n- k ak 


ak 


1 - 


a 2 n_k 


= f 


k 


2 n ~j-i 


- 1 


i - :i.f 


k 


2 n ~ 1 


= f (2 J+1 x — 1)[1 — jf (2x)\. 


The extension of g(x) to a continuous function on [0,1] follows from the facts that / extends 
to a continuous function on [0,1] and /( 2~i) = 1 /(j + 1). □ 
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Figure 3. Singular function associated with {6 n } 


A scaled version of the graph of y = f(x) appears in the graph of y = g(x). 

The restriction of g{x) to [1/2,1] is just a scaled version of f{x ): 

Corollary 3.10. g(x) = f(2x — 1) for x G [1/2,1], 

Recall that every positive rational number appears exactly once in the set {afc+i/a*; : k G N}. 
We prove an analogue for the sequence { 6 ^}. We use the expression “A = □” to mean that 
A = n 2 for some integer n. 

Theorem 3.11. Every element of {(a, 6 ) G N 2 : 4a6 + 1 = □} appears exactly once in the 
sequence {(frfc, 6 fc+i) : k G N}. 

Proof. Consider the following analogue of the (slow) Euclidean algorithm. 


M ffi : (a, b) 


( (a, a © b ) 
< (a © 6 , b ) 


^stop 


if a < b, 
if b < a, 
if a = b. 


Suppose (a, b) G N 2 , with 4a6 + 1 = □. If a © b < 0 then it is easy to see that (a — b) 2 < 1 
and thus a = b. In this case, since 4a 2 + 1 / □ unless a = 0, the only possibility is a = b = 0. 
Hence, M®((a, 6 )) G N 2 and, if this algorithm terminates at all, it must terminate at (0,0). 
With (a, b ) G N 2 , with 4a6 + 1 = □, let k := y/4ab + 1. If 0 < a < b, then a 2 < ak and thus 


a(a Qb) = a(a + b — k) = a 2 + ab — ak < ab. 


In general, the product of numbers in M^((a,b)) is strictly less than the product ab and 
thus the algorithm will eventually reach, without loss of generality, (0,6). If 6 = 0 then the 
algorithm stops. On the other hand, if 6 > 0, it is easy to see that M®((0, 6 )) = (0, 6 — 1), and 
thus the algorithm will terminate at ( 0 , 0 ). 

Let B n := ( 6 n , 6 n+ i). By the definition of the sequence { 6 ^}, it’s easy to see that for n > 1, 


A/® . B2m B2n+1 1 I B n 

and, moreover, if M® : (a, 6 ) i-G B n , then either (a, 6 ) = i? 2 n or (a, 6 ) = I? 2 n+i- 

If (a, 6 ) G N 2 , with 4a6 + 1 = □ is not of the form B n for some n, then all of its successors 
under M®, including (0,0), are not either - a contradiction. Hence every (a, 6 ) G N 2 , with 
4a6 + 1 = □ is of the form B n for some n. 

The pair (0, 0) appears only once and, in general, no pair appears more than once in {B n } 
for, otherwise, there exists a smallest n > 1 such that B n = B m for some m > n. Applying Af® 
to both B m and B n forces [ra/2j = |_m/2j and therefore m = n + 1. Thus b n = 6 n+ i = 6 n+ 2 , 
a contradiction. □ 
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A generalization of © is as follows: For a given number N, define 


x®y := x + y + JAxy + N. 

N 

Remark 3.12. a, b, a® b solve 

N 

2(x 2 + y 2 + z 2 ) - (x + y + z) 2 = N. 

Defining © in the obvious manner, 

N 

(a © b) © b = a. 

N N 

Every non-zero complex number z can be represented uniquely as re l6 for some positive r 

and some 9 € [0,27r) and so we define yfz := sjre ld I 2 . Hence © and © are well defined for 

N N 

complex N. 

We may then generalize {bk}- 

Definition 3.13. Given a (complex) number A, let c\ = C 2 = A and, for n > 1, 

^2 n 

^2n+l = Cn © Cn+l* 

N 

It turns out that such a sequence can be expressed as a linear combination of the sequences 
{a 2 } and {bk}- We first need a lemma. 

Lemma 3.14. For k > 1, 

a kbk+i + al+i&fc + 1 = Q-k a k+iy/ 46fcf>fc+i + 1- 

Proof. Let si, := \J Ab) ; bk+\ + 1. Note that 

b2k+i = bk + bk+\ + Sk- 


Then 


and so 


Similarly, 


and so 


Note that 


s 2k ~ + 1 — 4 bk(bk + frfc+l + Sk) + 1 

= ^b% + s\ + 4 bf.Sk = (2 bk + s k ) 2 

s 2k = 2 bk + Sk- 

s 2k+l = 4&2fc+1^2fc+2 + 1 = 46fc + i(6/ c + bk- 1-1 + Sk) + 1 
= 4 b 2 +1 + s\ + 46/ c+1 Sfc = (26fc + i + Sk ) 2 

S 2 fc+i = 26fc + i + Sfc. 

a \ b 2 + g -2^1 + 1 = 1 = a \ a 2 \/ 46162 + 1 


and so the lemma holds for k = 1 . 
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Suppose the lemma holds for a particular k. We show it works for 2k and 2k + 1 and thus, 
by induction, the lemma will be shown. 

a 2fe^2fe+l a 2fe+1^2fe + 1 = + bk+l + Sk) + (Ofe + afc+l) 2 6fc + 1 

= a 2 bk + a|6fc + i + a 2 Sk + a\bk + 2ak<ik+ibk + o, 2 +1 bk + 1 
= a\(2bk + s k ) + 2afcOfc_)_i6fc + a\bk +1 + a\ + ibk + 1 

= a|(26fc + Sk) + 2akCik+ibk + Qk a + k + lsfc 

— + Ofc+l)( 2 &fc + Sk) — 0 ' 2 k < ^ 2 k+ 1 ^ 2 k 

and thus the lemma works for 2k. 


a 2fc+1^2fc+2tJ2fc+2^2fc+l + 1 — (afc + afc+l) 2 ^fc+l + a k+l(bk + 1 + Sfc) + 1 

2 2 2 2 2 

= a kbk +1 + ^dkClk+lbk+1 + Clk+lbk+1 + + dk+lbk+l + Ofc+lSfc + 1 

= a\(2bk+i + Sk) + a\bk+i + a| + i bk + 1 + 2ak<ik+ibk+i 

— ®fc(“^fc+l “I - Sk) ~\~ OjkO'k+lSk “b 2 ofcCJfc+i&fc-|-i 
= ak+i{ak + ak+\){2bk+i + Sfc) = a2fc+202fc+i-S2fc+i 
and thus the lemma works for 2k + 1. □ 

Theorem 3.15. Given A,B, let Ck := Aa 2 + Bbk■ Then { Ck } has c\ = C 2 = A and, for 
N = 4AB + B 2 , 

C2n — 


Proof. 


and so 


and thus 


C2n+1 — C-n © Cn+1- 
N 


CfcC fc+ i + AH = (Aa 2 , + Bbk)(Aa% + i + + AB 

= A^aj(a 2 , +1 + B 2 bkbk+\ + AB(a\ +l bk + a\bk+\ + 1) 

= A^a|a 2 , +1 + B 2 bkbk +1 + ABakOk+w/ 4bkbk+\ + 1 

+ A = 4A 2 a 2 a^_)_ 1 + 4B 2 bkbk+\ + H 2 + 4ARafca; c _|_iy / 4bkbk+i + 1 
= (2^4a/ c a/ c _ ) _i + B\/ 4bkbk+i + l ) 2 

Cfc © Cfc+i = (Aa| + Bbk) + (Aa 2 , +1 + Bbk+i) + \/4cfcCfc + i + N 

= Aa 2 + Bbk + Aa\ +1 + Bbk +1 + 2AakOk+i + B\J 4bkbk+\ + 1 
= A(afc + Qfc+i ) 2 + -S(^fc + bk +1 + 4b i-b k +1 + 1 

= Aa^fc+i + H& 2 fc+1 = C2fc+1- 


Since 


c 2 fc — Aa| fe + Bb-jk — Aa\ + Bbk — c ki 


the theorem is shown. 


Example 3.16. Let N = —3, c\ = C 2 = 1, we see that A = 1, B = —1, and thus Ck = a? — bk. 
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If to every local cut point P in the fractal CP appearing in figure 2 one attaches a sphere 
above but tangent to the plane at that point with curvature ( 1 /radius) equal to the sum of 
the curvatures of the two circles meeting there, then one gets a 3-dimensional generalization 
of Ford circles. The curvatures (similarly, the product of local cut points and corresponding 
curvatures) along any circular arc are from a sequence {c n } for appropriately chosen N (see 
HB and references therein for a discussion of various types of “Ford spheres”). 

Consider the sequence {bk} written in tabular form: 

0 

0 1 

0 2 12 

0 3 2 6 1 6 2 3 

0 4 3 10 2 15 6 12 ... 


It is apparent that every column is an arithmetic sequence and, moreover, the defining differ¬ 
ences are respectively 

0,1,1,4,1,9,4,9,..., 

the squares of Stern’s diatomic sequence { aj .}. This is, in fact, true. We shall express this 
result as a formula. 

Theorem 3.17. For 0 < k < 2?, 

fyy'+ 1 +fc = a l + fyy'+fc- 

Proof. Assume 0 < k < 2 J . Since 2 J < 2 J + k < 2 J+1 , Theorem 3.6 implies 

bz>+ 1 +k = ^23+2 j+k = a 23+ 1 -{23 +k) a 23 +k-13 = a 2 3-k a k- 

By j9j formulas ( 2 ) and (3)] , 

CL2i+l-k = flfc + a 23-k 

and thus 

b 2 j+2 + k = a 23+ 1 -k a k = { a k + a ~ k ) a k = a \ + a 23-k a k = a k + ^2 3 +1 +k' 

The result follows by induction. □ 

Remark 3.18. {b 2 k- 1 } appears as [121 A119272], the product of numerators and denominators 
in the Stern-Brocot tree. 

Remark 3.19. For a fixed (x,y), z = x® y and z = x © y are the two solutions of 

2(x 2 + y 2 + z 2 ) - (x + y + z) 2 = 1. 


4. Fibonacci representations 

A Fibonacci representation of a number n is a way of writing that number as a sum of distinct 
Fibonacci numbers. One such representation is, of course, the Zeckendorf representation which 
is gotten by the greedy algorithm and which is characterized by having no two consecutive 
Fibonacci numbers. In general, a given n has several Fibonacci representations, the number 
of such we call R n . The sequence {R n } is extremely well studied; see papers by Klarner [7], 
Bicknell-Johnson m ei, and Stockmeyer B3I. for example. 
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A string of Os and Is is a finite word with alphabet {0,1} (equivalently, an element of 
{0,1}*). Often we denote such a word by u>. We shall think of such strings as Fibonacci 
representations: we shall assign a numerical value [w] to a string u by the formula 

= y ijFk+2-j- 

For example, [0100] = [0011] = 3 and [01010011] = 21 + 8 + 2 + 1 = 32. 

The generating function for {R n } has an obvious product formulation. 

Proposition 4.1. The sequence (R n ) satisfies 

oo oo 

y R n x n = (l +x Fi ) 

n =0 i =2 

where F n denotes the nth Fibonacci number. 


Next, we define the Fibonacci shift : 

p(n) := [ncf> + l/fi\ 

that satisfies p([w]) = [w0] for every string u. This shift has been studied before; for example, 
it appears in [6j graffiti, p. 301]. 


Theorem 4.2. For Ci E {0,1}, i = 2, 


TV 


TV 


p y = y Ci F i+1 . 

\»=2 ) i= 2 

Proof. By Binet’s formula (3), 

<\>F n = F n+ 1 -^ n . 

For any choice c* E {0, 1} for i = 2,..., IV, note that 


oo N 

r2n+l 


!/^ 2 =y y +1 < y^ < y y = 




i=2 


n=l 


and therefore 


Hence, 


TV 


o < — y cfif) — 4 > < i. 

i=2 


AT 


p y ctf; 


, i=2 


AT 


■y ^ - 

i=2 

y^ cj-^T+i + 


i=2 


TV 

- y cd? - 

. i=2 


N 

^ ^ Q-^i+1 • 


i=2 


□ 


In terms of p(n), we may define {!?„} recursively. Clearly, Rq = R\ = 1. A representation 
of n either ends in 0 in which case n = [cuO] where p{[oS\) = n or else it ends in 1 in which case 
n = [wl] and so n — 1 = [w0] = p([w]). Hence, for all n > 1, 

Rn := y Ri- 

p(i)e{n,n- 1} 
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Note that the function p 2 {n) := p{p(n )) = [ncjf 2 + l/^J is an example of a Beatty sequence 
(i.e., of the form [an + 6 J) and so has a complementary Beatty sequence, namely T(n ) := 
[rnji + 2/^>J . For example, 

p 2 (n) = 0, 3, 5, 8,11,13,16,18, 21, 24,... 

and 

T(n) = 1, 2 , 4, 6, 7,9, 10,12,14,15,... 

The following characterization could be used as a new definition of {P n }. 


Theorem 4.3. For n > 1. andT(n ) := [n(fr + 2/(j)\, 


and 


Rp2{n) Pn T Rn— 1 


-frr(n) — Rn- 


Proof. Since 4> £ (1,2), p(?r) € {p(n + 1) — l,p(n + 1) — 2}. Since 2cj> >3, (n — T)<f> + l/4> < 
(n + 1)0 + 1 /0 — 3 and so p(n — 1) < p(n + 1) — 2. Note that 

T[n) = [ncj) + 2/<j>\ = [{n + l)<f> + 1/0J - 1 = p(n + 1) - 1 

and therefore 

Rt{u) = ^ Ri = R n - 

p(i)e{p(n+l)-l,p(n+l)-2} 

We show the first equation in the theorem by a counting argument. By the definition of 
p(n), P 2 (n) = p{n) + n and so 

p 2 {n + 1) - p 2 (n) = p(n + 1) - p(n) + 1 G {2, 3}. 

For a given n, if n = [cj] then p(p(n)) = [wOO] and p(p(n + 1)) equals either [wlO] or [wll]. 
Suppose p 2 {n + 1) — p 2 (n) = 2 . The map ui i->- wOO is a bijection from representations of 
n to the representations of p 2 (n) ending in 00 while the map u i— > a ;10 is a bijection from 
representations of n — 1 to the representations of p 2 (n) not ending in 00. Hence the first 
equation holds. 

A similar argument holds when p 2 (n + 1) — p 2 (n ) = 3. □ 


Remark 4.4. The sequence {i? n } is thus analogous to the alternative form of Stern’s sequence: 


a 2n — a n , a 2n —i — a n T a n —±. 

For every word uj := uJoLJi...uj n € {0,1}*, we let |w| := n + 1 denote the length of u and 
define a point in the complex plane 

n 

P(u) :=^r fe (2w fc -l-i). 

k =0 

We form a graph G by putting an edge between P{oj ) and P(ooj) for j = 0,1, cu G {0,1}*. 
This graph is illustrated in Figure 4 below. Note further that P(oS) = P( uj') iff |w| = \u'\ 
and [w] = [u'\. Hence, we may consistently assign the integer [<u] to each vertex P(u) of the 
graph. This shows that is the number of downward paths from P(*) to P(w) and the 
graph can be thought of as a kind of hyperbolic Pascal’s triangle. In fact, the portion between 
0,01,010,0101,... and 1,10,101,1010,... is really just the “Fibonacci diatomic array” appearing 
in [2j. 

For v a vertex of the Fibonacci representation graph, let [v] be the number of downward 
paths from the top vertex to v. 
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Figure 4. Fibonacci Representation Graph with words in {0,1}*. 

Lemma 4.5. Along the nth row of the graph G, the function [u] forms an increasing sequence 
of consecutive integers 0,..., F n+ 2 — 2. 

Proof. Iterates of p(n) +1 starting at 0 yields the sequence 0,1, 3,6,11, F n+ 2 —2,... (provable 
by induction). Hence the last value of [u] in the nth row is F n+ 2 — 2. Since p(n + 1) — p(n) € 
{1,2}, the lemma follows. □ 

A consequence is the following surprising formula: 

p(p{p(n) + 1))) = p{p{p{n)) + 1) + 1. 

This graph has numerous interesting properties: 

• Every quadrilateral in the closure of the graph is either a square or a golden rectangle. 

• All the squares (actually hexagons) are congruent in hyperbolic space with area ln<^ 
(and, as hexagons, each edge has length In 0 ). The figure is thus an aperiodic tiling of 
part of the upper half-plane H (and can be extended to all of H ) where all the tiles 
are congruent! 

• The points along any row, when embedded in R form part of a one-dimensional qua¬ 
sicrystal. The lengths of the segments, appropriately scaled, form a word: 4>, 1, </>, 1, 

the “Fibonacci word”. 

• The vertices form a quasicrystal in H . 

• The graph is the Cayley graph of the “Fibonacci monoid” (a, b\abb = baa). 

• The graph can be constructed by the following recursive procedure starting with a 
single vertex; from each of the latest generation of vertices, draw two edges going 
southeast and southwest respectively, connect if a hexagon can be formed. Repeat. 

Something new with respect to the study of {R n } is the development of an analog of 
Conway’s box function. For k < F n _ 1 , define 

q(k, F n ) .— Rk/RF n +k. 

Lemma 4.6. For k = 0,..., F n _ 1 — 1, 

q(T(k),F n+ 1 ) = q(k,F n ) 

and 

q(p 2 (k), F n+2 ) = q(k, F n ) * q{k - 1, F n ) 
where * denotes “mediant addition”. 
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Proof. Note that 
and so, if k < F n ^\ — 1, 


T(n) = p(n + 1) — 1 
T{F n + k) = F n+1 +T{k). 


Then 


q(k,F n ) 


Rk _ _ 
RF n +k R 
Rr(k) 

R F n+ i+T(k) 


Rr{k) 

T(F n +k) 

= q(T(k),F n+1 ) 


and the first equation follows. Similarly, 


q(k,F n ) * q(k - 1 ,F n ) 


Rk Rk—l Rk F Rk —1 

- * - = - 

RF n +k RF n +k -1 RF n +k + RF n +k -1 
^P2(k) _ Rp2{k) 

Rp2(F„+k ) RF n+ 2+p2(k) 

q(p 2 (k),F n+2 ) 


and the second equation follows. 


□ 


As a consequence, if, as n —>• oo, k/F n converges to x G [0,1/0], then q(k, F n ) converges to 
some value, say Q(x). The function Q : [0, 1/0] —>• [0,1] is increasing and continuous. 



Figure 5. Analogue of Conway’s box function 

Note, however, it is not strictly increasing. 

Lemma 4.7. For j = 0, ...,F n - 1 — 1, 

RF n+ 2+j = RF n +j + Rj- 
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Theorem 4.8. The inverse of Q satisfies, on its irrational points of continuity, 


Q H x ) = X] 


(- 1 ) 


fc+1 




Cl+C2 + ...+Cfc)—1 


where x has continued fraction decomposition x = l/(ci + l/(c 2 + 1/(03 + ...))). 
Proof. Recall that Rp m +k = RF m+1 -k and so 

' ^'Pn +1 ^ _ ( jji IT? \ 

n + q(k.F m ) ~ R i + nR F ^ t ~ ~ 

Letting k/F m —>• x where x is a point of continuity of B , we see that 

1 


n + Q(x) 


= Q 


— x 

rk^Tl 


We may then rewrite: 


<t>- Q H®) _„-i{ 1 

^2n ^ 


n + x 


and the theorem follows. 


□ 


The function Q(x) extends past l/<f> but is no longer monotonic. 



Figure 6 . Analogue of Conway’s box function; larger domain 
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Patterns can be found by looking at the “crushed array” which is found by stacking rows 
of terms Rf^- i, Rf u+1 -2 sliding terms to the left on rows: 

1 

1 2 

1 2 2 

1 3 2 2 3 

1 3 3 2 4 2 3 3 

1433524425334 


The kth. column satisfies: x n +2 = x n + c with common difference c = 7?,/ ; ( Rq = 0). 
Alternatively, x n+ \ = x n + x n _i — x n -2 (a “dying rabbit” sequence). 


1 — T 1 %n—2 

Characteristic polynomial factors x 3 — x 2 — x + 1 = (x — l) 2 (x + 1) so every example is of the 
form x n = a + bn + c(—l) n . Hence, {x 2 n} and {x 2 n+i} are arithmetic sequences. 

Xn +1 — T X n _i Xjj —3 

e.g., [12^ A023434] x 4 — x 3 — x 2 + 1 = (x — l)(x 3 — x — 1), so every example is of the form 
a + 6r” + ct 2 + dr^ where r\ is the “plastic constant”, 1.32471795..., the smallest Pisot number, 
and r 2 ,r 3 are its algebraic conjugates. Such examples are always a constant plus a Padovan 
sequence y n+ \ = y n _i + y n _ 2 . E.g., [HI A000931] 


Xn+l — “h 3'n —1 


n—1) 


is always a constant sequence. 


5. Extending Binet’s formula 

Let sjr(n) be the number of terms in the Zeckendorf representation of n (e.g., sf{ 27) = 3). 
Equivalently, s_p(n) is the least number of Fibonacci numbers that sum to n. This sequence, 
for n = 0,1,..., is [A007895] and starts 


0,1,1,1,2,1,2, 2,1,2,2,2,3,1,2, 2, 2,3, 2,3,3, ... 

Using notation of the previous section, we see that sjr(n) satisfies the recursion: 

sf([w0]) = s f ([w]),sf([w01]) = sf(M) + 1 


which translates to 


s F (p(n)) = s F (n), s F (p 2 (n) + 1) = s F (n) + 1 
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where p(n ) is the “Fibonacci shift” defined in Section 4 (just after Proposition 4.1). The 
crushed array for this sequence is 

1 

1 

1 2 

1 2 2 

1 2 2 2 3 

1 2 2 2 3 2 3 3 

1222323323334 


Note that columns are constant and that the limiting row is s F (n) + 1. 

Replacing s 2 (ro) by s F (n) in Binet’s formula for Stern’s sequence yields our third variant of 
Stern’s sequence: 

n 


c n+1 

k =0 


The sequence starts, for n = 1, 2,..., 


1,1, 2, 3, 2,4,3, 3, 6 ,4, 6 , 6 ,4, 8 , 6 , 7,10, 6 , 9, 7,5,11, 8 ,.... 

It is always integral since c n +1 is an algebraic integer in Z[a] invariant under complex conju¬ 
gation. 

A crushed array for this sequence is: 

1 

1 

2 3 

2 4 3 

3 6 4 6 6 

4867 10 697 

5 11 8 11 13 8 14 10 9 15 9 13 11 

7 15 11 15 19 12 19 14 11 21 14 19 19 

The first column, x n := { c Fn } apparently satisfies the Padovan recurrence: x n+ 2 = x n + x n -\. 
Moreover, every column is apparently a “dying rabbit” sequence: x n+ \ = x n + x n -\ — x n s 
or, more precisely, if x n := c Fn+k + c k , then x n +2 = x n + x n -\. This is indeed the case which 
we now prove. 


Theorem 5.1. For k < F n _ 2 , c Fn+2+k = c Fn+k + c k + c Fn _ 1+k . 

Proof. Given a string X of integers, let X denote the reverse of string X, let X + denote the 
string X with 1 added to every integer, and let X~ denote the string X with 1 subtracted 
from every integer (e.g., if X = 1223, then X = 3221, X + = 2334, and X = 2110). If 
X := let G{X) := o atj - Finally, given strings X , Y. we let XY denote the 

concatenation of the two strings and X — Y denote the pointwise difference (e.g., if X = 457 
and Y = 123 then XY = 457123 and X — Y = 334). 

Let s n := s F (n) be the number of terms in the Zeckendorf representation of n. For any 
interval /, let s/ denote the string Si 1 Si 2 ...Si k where i\ < i 2 < ... < i k and {ii,i 2 , ...,*&} = 7nN. 
Then c n = G((s[ 0jn ) — S[o,n))) where the difference of two strings is the string of differences. 
Since we will use this formula, we let A/ = sj — sj so that c n = G(A[ 0 r) j). 
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By the definition of s, it’s clear that S[_F ni _F„+fc) = s [o k ) ^ ^ — F n _ i- Since 
s [0,F n +k) = s [0,k) s [k,F„) s [F„,F n +k) = s [0,k) s [k,F n ) s ^ k y 


it follows that 

and thus 

Hence, because a + a = 1, 


s [0,F„+fc) S {o,ik) S [fe,-Fn) S [0,fe)> 
A [0,F„+fc) = A [0 ,fc) A [fc,-Pn) A J,fc)- 


CF„+fc = O’ 1 Cfe + G(A[ fc;Fn )) + CTCfe = Cfe + G(A[ fcii r n )). (*) 

Assuming k < F n _ 2 , we see that 

s [0,F n +2+fc) s [0,k) s [k,F n ) s [F n ,F n +k) s [F n +k,F n +i) s [F n -\-i,F n -\-i+k) s [F n -\-i+k,F n +2) s [F„+2,F n +2+k) 

and thus 

«[0,F n+2 +fc) = ^io,fe)^ife,F n )^jo,fe)^ife,F n _i)^[0,jk)®[fe,^n)^[0,fc)- 

Hence, 

A [0,F„ +2 +fc) ~ A [0,fc) A [fc,F n ) A [0,fc) A [fc,FVi— 1 ) A [0,fc) A [fc,FVi+A:) A [0,fc) ‘ 

Applying G: 

CF n+2 +k = & 1 Ck + cr 1 G(A[ fcii 7 n )) + Cfc + G{ A[ fcjFn _ 1 )) + c k + crG(A[ fc!Fn )) + ac k - 
Again, since cr -1 + a = 1, and by (*), we have 

cf„ +2 +A: = 3 c k + G(A[ kiFn )) + G(A [fc)F „_ 1 )) = c k + c Fn+k + c Fn _ 1+k . 

□ 


There are many patterns in the crushed array. Two such patterns can be proven by induction 
based on the previous theorem. 

Corollary 5 . 2 . c Fn + c Fn _ 1+ 2 = c Fn+ 1 and c Fn+ 1 = CF n+1 +2 /or all n. 

We have many other questions or apparent properties, all waiting for a proof (though, of 
course, of varying difficulty). 

• Five inequalities: c <72 ( n)+ 1 > c^j > c^j > c a(n) >c n > 0. 

• The minimum of each row in the crushed array is the leftmost element. (If true, then 
the last inequality above, c n > 0, is true). 

• If Cn A 0 for all n, then what do these numbers count? 

• The following sequences have crushed arrays with columns satisfying x nF \ = x n + 
x n -i - x n -j for given j: 

{s F (7i)} has j = 1, 

{R n } has j = 2, 

{Cn} has j = 3. 

Is there a general principle at work in this progression? Is there a similarly defined 
sequence with j = 4 for example? 
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